Accurate Group Delay Measurement for Radial Velocity Instruments Using the 
Dispersed Fixed Delay Interferometer Method 
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ABSTRACT 

The Dispersed Fixed Delay Interferometer (DFDI) method is attractive for its low 
cost, compact size and multi-object capability in precision radial velocity (RV) mea- 
surements. The phase shift of fringes of stellar absorption lines is measured and then 
converted to an RV shift via an important parameter, phase-to- velocity scale (PV scale) 
determined by the group delay (GD) of a fixed delay interferometer. Two methods of 
GD measurement using a DFDI Doppler instrument are presented in this paper: 1), 
GD measurement using white light combs generated by the fixed delay interferometer; 
2), GD calibration using an RV reference star. These two methods provides adequate 
precision of GD measurement and calibration given current RV precision achieved by 
a DFDI Doppler instrument. They can potentially be used to measure GD of an fixed 
delay interferometer for sub-meter precision Doppler measurement with a DFDI instru- 
ment. Advantages and limitations of each method are discussed in details. The two 
methods can serve as standard procedures of PV scale calibration for DFDI instruments 
and cross checks for each other. 

Subject headings: methods: planetary systems-techniques: radial velocities 
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1. Introduction 



As of Apr 2012, there are over 700 discovered exoplanets, and most of them are d etected by the 
radial velocity (RV) techn iqugj- RV precision of 1 m-s~^ has been routinely ach ieved ( Bouchv et al. 



2009 



Howard et aPboiol ) with instruments such as HARPS (JMayor et ahlbood ) and HIRES (|Vogt et al 



19941 ). which are cross-dispersed echelle spectrographs. While cross-dispersed echelle spectrographs 
are commonly used in instruments for precision RV measurements, a method using a dispersed 



fixed delay interferoni eter (DFDI) has offered an alternative method (jGe et al. 



2010 



Lee et al. 



2011 



2006 



Fleming et al. 



In this method, a Michelson-type interferometer is used in combination 



^http://exoplanet.eu/; http://exoplanets.org/ 



with a moderate resolution spectrograph, RV signals are then extracted from phase sh i ft of in- 



terfer ence fringes of stellar absorption lines (JErskine &: Gd l200d : iGe et al. 



20031) ■ The deta i ls abo ut the DFDI theory and applications are discussed in 



2002; Ge 



2002 



Ersking 



van Evken et al.l (J2010l ) 



and I Wang et al.1 (|201ll ). Instrument adopting the DF DI method has denionstrated advantages such 
as low cost, compact size and multi-obj ect capability (JGell2002l : iGe et al.ll2006l : [Fleming et al.ll20ld : 



Lee et al. 



2011 



Wisniewski et alJl2012l ). 



In the DFDI method, a fixed delay interferometer (jWan et alj l2009l . l201lll plays a crucial 
role in creating stellar spectral fringes for high precision RV measurements (JGd |2002| : lErskine 
20031 ). The Do ppler sensitivity ca n be optimized by carefully choosing the group delay (GD) of the 
interferometer (jWang et al.ll201ll ) . More specifically, GD of an interferometer should be chosen such 
that the spatial frequency of white light combs (WLCs) matches with that of a stellar spectrum 
after rotational broadening. GD is defined by the following equation: 



GB{iy) 



1 

'2^ 



di/' 



(1) 



where 4> is phase shift and u is optical frequency. The interferometer in a DFDI instrume nt is usually 



designed to be field-com pensated to minimize the influence of input beam instability ()Wan et al. 



2009 



Wang et al.ll2010l ). It is realized by carefully selecting glass materials and thicknesses of two 



second surface mirrors such that their virtual images are overlapped. Because glasses are used in 
the optical paths, (j) does no longer linearly change with frequency, therefore GD is dependent on 
optical fr equency. An inaccurate GD measurement may significantly limit the RV measurement 
accuracy (JBarker fc Schulerlll974l ; Ivan Evken et al.ll2010l ). 



In practice, there may be several methods of measuring GD: 

1. Calculate GD based on glass refractive indices using Sellmeier equation and thicknesses from 
manufacturer specification. 

2. Forward model the spectrum of a known spectral source, such as an Iodine cell or a Th-Ar 
lamp. 

3. Measure phase and frequency using a while light source, such as a tungsten lamp. 

4. Calibrate GD using a source with known velocity. 



Method 1 is straightforward but may lack of adequate precision because of uncertainty in 
parameters in Sellmeier equation and manufacturer tolerance for glass thickness. Method 2 holds 
great promise for accurately determining GD but there some current practical issues preventing us 
from adopting this method (see more detailed discussion in §5.2.ip . We will use Method 3 and 4 
to measure GD of an interferometer in this paper. 

In the DFDI method, GD determines the phase-to-velocity (PV) scale, the proportionality 
between the measured phase shift and the velocity shift. Since the DFDI method is realized by 



coupling a fixed delay interferometer with a post-disperser, the resulting fringing spectrum— stellar 
absorption lines superimposing on the WLCs— is recorded on a CCD detector (illustrated in Fig. 
[1]). The fringe phase is expressed by the following equation: 

H'^,y) = , (2) 

c 

where y is the coordinate along slit direction, which is transverse to dispersion direction, r is the 
optical path difference (OPD) of an interferometer and c is the speed of light. Two mirrors (arms) 
of the interferometer are designed to be tilted towards each other along the slit direction such that 
several fringes are formed along each z^ channel. The intersection of a stellar absorption line and 
a WLC moves (from Pq to P in Fig. [1]) if there is a shift of an absorption line due to a change of 
stellar RV. Consequently, a small change of (j) in the dispersion direction, A(j)x, is induced: 

A , d4> X d(f) dv . 

dv dv dv 

d(j) v . „ . 
= --•-• Aw = r • Av, (3 

dv c 

where F is defined as phase-to- velocity scale (PV scale). It is determined by the GD of an interfer- 
ometer, which becomes explicit if Equation ([T]) and ([3]) are combined: 

F = -27r-GD--. (4) 

c 

At resolutions typically adopted by the DFDI method (5,000 < R < 20,000), stehar hues 
(line width~0.lA) are not resolved and a measurement of A(j)x is extremely difficult. Instead, Acpy, 
phase shift along y direction can be measured, which is equal to Acj)^ if the combs generated by 
an interferometer are parallel to each other. This is a good approximation at very high orders of 
interference. The advantage of measuring A(f)y instead of A(px is seen from Fig. [H in which the 
physical shift in the u direction is amplified in y direction, the amplification rate is determined by 
the relative angle between the interferometer combs and a stellar absorption line. Therefore, A(py is 
relatively easier to measure compared to (px and it is measured by fitting a well-sampled periodical 
fiux signal along the y direction in the DFDI method. Compared to conventional high-resolution 
Echelle method, the number of freedom for the DFDI method in the fitting process is much less 
and small Doppler phase shift can be relatively easier detected with a simple functional form, i.e., 
a sinusoidal function. However, we want to point out that while the DFDI method provides a 
boost in instrument Doppler sensitivity, the Doppler sensitivity is not strongly dependent on the 
amplification rate because fiux slope decreases as amplification rate increases, which negates the 
gain of phase slope. 

The paper is organized as follows: in ^ we present a new method of GD measurement using 
a DFDI Doppler instrument, which exploits WLCs generated by a fixed delay interferometer. In 
^ we present another method of GD calibration using an RV reference star. Observation results 
after implementing the newly measured GD are presented in ^ In f|5l we summarize and discuss 
the new results in this paper. 



2. GD Measurement Using White Light Combs 
2.1. Method 

MARVELS (Multi-object Apache Point Observatory Ra dial Velocity Exoplanet Large-area 



Surv ey) is part of the Sloan Digital Sky Survey (SDSS) III (JGunn et al.l l2006l : lEisenstein et al 



I2OIII ) . The instrument covers a wavelength range from 500 nm to 570 nm and us es a post-disper sive 
grating with a spectral resolution of 11,000 after a fixed delay interferometer (JGe et al.ll2009l ). A 
Th-Ar emission lamp and an io dine absorption cell serve as wa velength calibration sources. The 



instrument setup of MARVE LS (IGe et al. 



measure GD as described in Kovacs et al 



I 



2009l:IWan et al.ll2009l') is similar to th e equipments that 



(|l995l ) and lAmotchkina et al.l (|2009l ) , in which a white 



light interferometer ( WLI) is combined with a post-disperser. However, the OPD is scanned by a 
moving picomotor in lAmotchkina et al.l (J2009l ) while it is realized by two relatively tilted arms in 
the WLC method using MARVELS instrument. WLCs are generated by the interferometer when 
fed by a white light source (e.g., a tungsten lamp). (j){i'), the phase of each frequency channel ly, 
is measured and then unwrapped to remove ambiguity of 27r. GD is then derived by taking the 
derivative of 0(i^) according to Equation ([T|). 



2.2. Data Reduction 

Standard spectroscopy reduction procedures are performed with an IDL data reduction pipeline 
dedicated to MARVELS. Fig [2] shows an example of normalized flux as a function of frequency for 
a processed spectrum. A zoom-in sub-plot shows the WLCs produced by frequency modulation 
of the interferometer. Visibility, defined as the ratio of half of peak-valley value to the DC offset, 
increases with frequency in the red part of the spectrum. The increasing visibility in the blue end 
of the spectrum is not physical but caused by an increasing photon noise and our algorithm of 
visibility calculation. 



^The fringe phases as a function of i^ are calculated by the Hilbert transform technique (JRochford &: Dyei 

I999I ): the signal H{i') is obtained by firstly removing the negative Fourier components of -F(i^)— the 



flux distribution with frequency— and then conducting an inverse Fourier transform. The phases 
(f){i') are obtained by calculating and unwrapping the arguments of H{v). We find that the phase 
change between pixels exceeds vr in the blue part and therefore the phase unwrap cannot be suc- 
cessfully applied, so we decide to use only part of the spectrum with a pixel range from 1800 to 
3800 for phase unwrapping. A third-order polynomial is used to fit (^ as a function v. GD is then 
calculated according to its definition (Equation ([T])). 



2.3. GD Measurement Results 

The top view and side view of the MARVELS interferometer are shown in Fig. [3l 60 fibers are 
mounted and each creates two spectra, one is picked from the forwarding beam and the other one is 
from the retm'ning beam (see Fig. [3|). In total, 120 spectra are formed, allowing us to measure GD 
at 60 positions on the interferometer along vertical (slit) direction. Each position corresponds to a 
fiber number. There are 24 pixels along the slit direction for each spectrum. We chose 15 rows in 
the middle to measure GD because of relatively higher photon flux, and thus smaller photon noise 
in the middle region of the spectrum. The top panel of Fig. S] shows phase measurement results 
for center row as a function of frequency at different fiber numbers. Phase fitting residual (shown 
on the bottom panel of Fig. [U RMS=0.9 rad) is consistent with photon-noise limited measurement 
error (see i l2.4l for details). GD for a particular fiber number is obtained by averaging the results of 
GD measurements for those rows associated with the fiber. Fig. [5] shows the results at i^=550 THz 
as a function of fiber number. Note that the two arms of the interferometer are intentionally tilted 
to each other and the 60 fibers are evenly mounted along the slit direction. The measured GDs 
should gradually vary with fiber number. We use a second-order polynomial to fit the GD variation 
with the fiber number. The fitting residual has an RMS of 0.0046 ps. Fig. [6] shows fitted GD as a 
function of frequency for different fibers. GD varies 0.15 ps (0.6%) across measurement range from 
540 to 565 THz. Ignoring GD dependence of frequency would result in 180 m • s~^ measurement 
offset between two ends of measurement range (assuming a true RV of 30,000 m • s~^, which is a 
typical stellar RV value due to the Earth's barycentric motion). Table [1] provides the polynomial 
fitting coefficients of GD vs. fiber number at different frequencies within measurement range. 



2.4. GD Measurement Error Analysis 

Two physical parameters, (j) &iid u, are measured in the experiment. The uncertainty of the 
(f) measurement is ~0.8 rad under photon-noise limited condition assuming a S/N of 120 and a 
typical fringe visibility of 1.5%. The uncertainty due to the wavelength calibration is ~0.002 
THz (O.O2A). We conduct a bootstrapping process to investigate the uncertainty of GD caused 
by the measurement uncertainties of (/> and v. We add gaussian noises with standard deviation of 
measurement errors to both (j) and v and calculate the group delay. We run 1000 iterations for 
bootstrapping in order to estimate the uncertainty of GD. The median of the relative error of GD 
measurements, 5GD/GD, is 4.4 x 10~^. In comparison, the median of the relative GD measurement 
error is ~ 1.8 x 10~^ after smoothing by fitting a polynomial to GD variation with the fiber number. 
This number does not agree with the relative GD error predicted by the bootstrapping experiment. 
We suspect that the uncounted error in the bootstrapping simulation comes from image distortion 
due to optics which the data pipeline has not fully corrected for, e.g., spectrum curvature, spectral 
line slant, etc. In a 2-D spectrum as illustrated in Fig. [H the phase shift between adjacent 
pixels along slit direction is ~0.6 rad, and the phase shift between each wavelength chanel is '~2.5 
rad. An imperfect spectrum curvature tracing tends to shift pixel in the slit direction while an 



imperfect slant correction can affect pixel shifting in both slit and dispersion directions. The range 
of unwrapped phase is ~4000 rad. For one fiber, if a gradually-changing phase error is introduced 
by the data pipeline when correcting for the optical distortion, for example, 0.4 rad deviation from 
true value at one end while no deviation at the other end, a relative error of GD would be caused 
with an estimation of 0.4/4000 = 1 x 10~^. If different fibers are treated independently, which 
is the case for the MARVELS data reduction pipeline, then this gradually-changing phase error, 
introduced by imperfect optical distortion correction, may explain the standard deviation error we 
see after the polynomial fitting for GDs as a function of fiber number. 



3. GD Calibration: Observing an RV Reference Star 

3.1. Method 

A deviated PV scale would result in an inaccurate velocity measurement given the same amount 
of fringe phase shift: 

/^(j) = Av' -r' = AvT, (5) 

where Av' represents a measured velocity shift while Av represents a true velocity shift. Combining 
Equation @ and ([5]), we obtain the following equation from which GD can be calculated by using 
the measured velocity shift of an object with a known velocity. 



This approach is similar to that of [Barker &: Schuleij (|1974| ). but the difference is that the latter 



applied correction for discrete laser frequencies while we seek corrections for a continuous frequency 
distribution. 

In order to realize the method of GD calibration using an RV reference star, we need to: 1), 
assume a GD' that is close to the true value of GD; 2), measure velocity shift Av' based on an 
assumed GD'; 3), know the true value of the velocity shift of an RV reference star. 



3.2. GD Calibration Precision 

The calibration precision using an RV reference star is determined by measurement error of 

Av': 

GD'-^v' 
6GD = -^^, (7) 



where 6v' is RV measurement uncertainty. I Wang et al.l (|201ll ) provided a method of calculating 



photon-limited RV uncertainty for the DFDI method. Under photon-noise limited condition, we 
expect the GD calibration error to be determined by the photon-limited RV uncertainty within the 
instrument band width, which is provided in Table [2j In Equation ([7|), GD' is usually estimated to 



be within a few percent of true GD, Av is statistically ~10,000 m • s~^ given a uniform reference 
star distribution and a quarter year observational availability. Therefore, relative error of GD 
measurement is -^ 2 x 10"'^ for an RV reference star with a T^s of 4500 K. 



4. Implementation of Measured GD in Astronomical Observations 

We use an RV reference star, HIP 14810 (V=8.5), as an example to show the RV measurement 
results after implementation of the newly measured GD using the WLC met hod. HIP 14810 is a 



star known to harbor 3 planets and its RV jitter is estimated to be 2 m • s (jWright et al.ll2009l ). 
After RV changes due to instrument drift, the Earth's barycentric motion and orbiting planets 
are removed, RV RMS error is 17.13 m • s~^ but has not reached the predicted photon-limited RV 
uncertainty (4.8 m • s~^, S/N=80 with a half wavelength coverage from 535 to 565 nm). RV RMS 
error is expected to be further reduced after the data pipeline is improved in the future. 

We also examine the reference star GD calibration method. We use one spectral block within 
measurement range centering at 550 THz (540-560 THz) and set GD' to be an arbitrary value 
of -23.873 ps. The measured RVs (barycentric velocity not corrected) are shown in Fig. [71 Af- 
ter applying correction according to Equation ([6]), we find that the GD is -25.107±0.027 ps. In 
comparison, GD measurement result of fiber number 51 (the fiber for HIP 14810) using the WLC 
method gives -25.091ib0.005 ps (refer to Table [1]). We confirm that the GDs measured by these two 
methods are consistent with each other at 68% significance level. 



5. Summaries and Discussions 

5.1. Summaries 

The PV scale is an important parameter in the DFDI method that translates a measured phase 
shift to an RV shift, and is determined by the group delay (GD) of an interferometer. We have 
provided and discussed two methods of GD measurement and calibration: 1), GD measurement 
using white light combs (WLCs) generated by the interferometer in a DFDI Doppler instrument; 
2), GD calibration using an RV reference star (RS). Table [3] summarizes the main results and the 
comparison between these two methods. The accuracy of GD meas urement is sufficient f or current 



RV precision achieved with instruments using the DFDI method (JFleming et al.ll2010l : iLee et al. 



2011 



Muirhead et al.ll2011l ). However, higher measurement and calibration precision is required in 
the near future as higher RV precision is achieved by DFDI instruments in search for exoplanets. 
RS and WLC methods can serve as complementary methods of GD measurement and calibration 
for DFDI instruments. 



5.2. Discussions 

5.2.1. White Light Corah (WLC) Method 

The GD measurement using WLCs created by the interferometer provides a direct way of 
calibrating the PV scale. In the region where combs are visible, effective S/N is relatively low 
(~15) because of low comb visibility (1.5%). In addition, GD cannot be measured in the region 
where combs are not visible, which limits the application of this method. We are able to measure GD 
in a region that accounts for half of the spectrum coverage. Extrapolation beyond the measurement 
range may result in large uncertainties. The major issue facing the method is that the data reduction 
pipeline may have introduced unknown errors while correcting optical distortions such as spectrum 
curvature and slant. 

In principle, we can use a tungsten lamp with an iodine cell or a Th-Ar lamp instead of a 
tungsten lamp in order to increase the fringe visibility. However, there are some practical concerns 
that hinder us from applying the above solutions: 1), line blending, because of low spectral resolu- 
tion, many spectral lines can not be resolved and it is not certain at this stage how line blending 
affects phase measurement; 2), illumination correction, which is required to correct for illumination 
profile in slit direction in order to properly measure the phase. We adopt a self-illumination correc- 
tion procedure in the pipeline which requires a certain continuum level to be successfully achieved. 
Th-Ar is less affected by line blending if a careful line selection process is involved, but it does not 
have enough continuum level for self-illumination correction. An experiment is being conducted 
in which a second interferometer is used to improve the visibility of WLCs so that GD is more 
precisely measured at a higher effective S/N for a wider frequency coverage. 

When compared to previous work in the field of GD measurement, lAmotchkina et al.l (120091 ) 



achieved a measurement precision of 1 x 10~^ ps. Our measurement of GD has a typical accuracy 
of ~ 6 X 10~^ ps (limited by effective S/N and systematic errors), which is more than an order 
of magnitude lower. However, it is shown in ^2.41 that substantial improvement would be able to 



be ach ieved once the data reduction pipeline has a better handle of optical distortion. I Wan et al. 



(|2010l ) measured GD for the MARVELS interferometer using a scanning WLI method and achieved 
a precision of 0.6 x 10~^ ~ 2.4 x 10~^ ps, which is more than two order of magnitude better than 
the results in this paper. However, there are practical concerns using GD measurement results 
from the scanning WLI method because they are not measured in situ, therefore it is not easy to 
associate a position in a WLI measurement to a fiber position. 



5.2.2. Reference Star (RS) Method 

The GD calibration using an RV reference star (RS) is a self-calibrating process and has the 
potential of achieving a high calibration precision if the following requirements are met: 1), the 
RV reference star has a large velocity shift during a observation window; 2), the RV reference star 



is bright; 3), the data reduction pipehne is able to produce the photon- hmited RV precision. In 
addition, GD is practically measured within a certain band width: 

GD(z.) = ^^^ \\y , (8) 

where u;(i^) is weight function. The band width, Av, should be small such that the dispersion effect 
is negligible. The limitations stated above prevent us from precisely determining the PV scale at 
the position of each fiber because not every fiber has a continuous observation on a bright known 
RV reference star. For MARVELS, the brightest RV reference star available has V mag of 8 and 
the resulting S/N is ~100 per pixel. However, the RS method is a very promising approach for a 
single-object DFDI instrument because only one bright reference star is required in the field. We 
are planning to apply this method in calibrating GD of another DFDI instrument (EXPERT) at 



KPNO 2.1m telescope (JGe et al.l 120101 ). Note that the S/N can be further increased by conducting 



multiple independent measurements and increasing instrument throughput. 
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Fig. 1. — Illustration of the DFDI method. Tilted lines represent interference combs generated by 
an interferometer. Vertical line represents an stellar absorption line (solid: original position with a 
frequency of fo! dashed: shifted position with a frequency of z^). 
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Fig. 2. — The normalized flux and visibility (7) as a function of frequency of a tungsten spectrum 
taken with MARVELS. The solid line is the normalized flux and filled circles represent visibilities 
in different frequency channels. 
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Fig. 3. — Top: top view of an individual fiber beam feeding of the MARVELS interferometer. Two 
spectra are formed by one fiber. One (Slit A) is from the returning beam arm while the other one 
(Slit B) is from the forwarding beam arm. Bottom: side view of the fiber array beam feeding of the 
MARVELS interferometer. There are 60 fibers yielding 120 spectra. Note the exaggerated wedge 
angle of the shown second surface mirror, GD gradually changes along the vertical direction. 
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Fig. 4. — Top: white light combs phase as a function of frequency at different fiber locations. 
Bottom: phase residual after third-order polynomial fitting. 
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Fig. 5. — Measured group delay as a function fiber number. Filled circles are measured results, 
solid line represents the best second-order polynomial fitting with an RMS fitting error of 0.0046 
ps. Measurement results can be found in Tabled] at other frequencies. 
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Fig. 6. — GD as a function of frequency at different fiber numbers. GD measurement results vs. 
frequency and fiber numbers can be found in Table [H 
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Fig. 7. — Top: measured {Av') and true (Af) RVs of HIP 14810 (barycentric velocity not corrected) 
over a period of 70 days. Bottom: the ratio of Av' and Av as a function of time. 
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Table 1. GD measurement results as a function of spectrum number 

{GD{#) = Co + C*! • # + C2 • #2) and standard deviation (5GD) at different 

frequencies (z^) 



zy[THz] 



Co 



Ci 



Co 



5GD[vs\ 



540.0000 
542.0000 
544.0000 
546.0000 
548.0000 
550.0000 
552.0000 
554.0000 
556.0000 
558.0000 
560.0000 



-2.5195975233e+01 
-2.5216854974e+01 
-2.5219400268e+01 
-2.5232735338e+01 
-2.5246276599e+01 
-2.5259435353e+01 
-2.5270775315e+01 
-2.5286851225e+01 
-2.5295872014e+01 
-2.5303654859e+01 
-2.5319347252e+01 



4.1224184409e-03 
4.8064043246e-03 
3.8524896642e-03 
3.5346857088e-03 
3.5171101494e-03 
3.3877963986e-03 
2.8029235840e-03 
3.5335302125e-03 
3.4000221425e-03 
2.3432563708e-03 
2.6464319719e-03 



-1.3537701123e-05 0.0066 

-2.3438315106e-05 0.0070 

-8.1355543995e-06 0.0057 

-2.6247780042e-06 0.0055 

-3.3394045215e-06 0.0047 

-1.7477210448e-06 0.0046 

9.6338400722e-06 0.0066 

-4.0412480867e-06 0.0093 

-2.5813725398e-07 0.0067 

1.8389109388e-05 0.0087 

1.2493738694e-05 0.0075 



Table 2: MARVELS predicted RV uncertainty (at an average S/N of 100) vs. T^^ 



%s[K] 


4500 


5000 


5500 


6000 


6500 


7000 


fo'[m-s-^] 


1.9 


2.3 
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3.5 


4.0 



Table 3: Comparison between two methods of GD measurement and calibration 







WLC 




RS 


Spectrum Coverage 




Half 




Full 


S/N 




-15 




100 for V~8'' 


Current precision 


4.6 


X 10-3 ~ 9.3 X 10" 


"3 ps 


0.027 ps 


Current RV error'' 




~ 2 m-s-i 




10.8 m-s-i 


Potential precision 




~ 3.1 X 10"'' ps 




--^ 5 X 10-3 ps 


Potential RV error^ 




~ 0.1 m-s"^ 




~ 2 m • s"^ 


Dependence on observation 




X 




/ 


Dependence on pipeline 




/ 




/ 













Note. — a: assuming MARVELS throughput; b: 
10,000 m • s"^ according to Equation (|6]). 



RV error is calculated assuming a true velocity shift (Au) of 



